Gauge Theory for the Rate Equations: Electrodynamics on a Network 
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Systems of coupled rate equations are ubiquitous in many areas of science, for example in the description 
of electronic transport through quantum dots and molecules. They can be understood as a continuity equation 
expressing the conservation of probability. It is shown that this conservation law can be implemented by con- 
structing a gauge theory akin to classical electrodynamics on the network of possible states described by the 
rate equations. The properties of this gauge theory are analyzed. It turns out that the network is maximally 
connected with respect to the electromagnetic fields even if the allowed transitions form a sparse network. It is 
found that the numbers of degrees of freedom of the electric and magnetic fields are equal. The results shed light 
on the structure of classical abelian gauge theory beyond the particular motivation in terms of rate equations. 
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Introduction. — ^Let us consider a system that assumes states 
\i) with probabilities Pi, where the states form a finite 
or countable set. The probabilities must add up to unity, 
J2i Pi — 1- If the rates of change of the probabilities are 
linear functions of the probabilities, we can write 
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(1) 



This is a set of rate equations, containing the transition rates 
Rij from state |j) to state The first term under the sum de- 
scribes transitions from other states to whereas the second 
describes transitons out of state Often, many of the Rij are 
zero. Importantly, Eq. (|T|) conserves probability: Pi — 0. 

Rate equations are ubiquitous in science, in particular in 
physics, chemistry, and biology. They describe systems far 
from equilibrium, such as lasers |1], semiconductor devices 
[2], quantum dots [3], chemical reactions [4], enzyme kinetics 
and biological populations ll^l- 
As a specific example, we discuss a molecular quantum dot 
coupled to leads. The time evolution of the complete system 
of molecule and leads is described by the von Neumann equa- 
tion for the full density operator, p = ~i[H,p\, where H is 
the Hamiltonian of the complete system and h — 1. If one is 
interested in the properties of the molecule, it is advantageous 
to integrate out the states of the leads to obtain an equation 
of motion (master equation) for the reduced density matrix 
p'^ in the many-particle Hilbert space of the molecule alone 



p"- in tiie many-particle 
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16l llTl uM- This requires 



approximations. For example, if the tunneling amplitude be- 
tween molecule and leads is small compared to the typical 
molecular energy-level spacin g, o ne can apply perturbation 
theory [1101 [m, 15, lelo 



1811 . The master equation 



contains diagonal and off-diagonal components of p'^. If one 
assumes that the off-diagonal components, i.e., the superpo- 
sitions, decay rapidly on the relevant timescale [19], one is 
left with coupled rate equations ([T]l for the diagonal compo- 
nents, which are just the probabilities Pi = pf^ of molecular 
many-body states. 

Our discussion is not specific to a particular incarnation of 
the rate equations. Rather, it starts from the general princi- 



ple of conservation of probability. We consider the network 
formed by the states The transitions between them are the 
edges of this network. The rate equations ([T]l can be written 
as a continuity equation on the network 1,20.1 . 



Pi — Jij , 



(2) 



where Jij is a probability current defined as an antisymmetric 
quantity on the edges Bill . 



T = R P — R P 



(3) 



The conservation law (|2]l suggests to look for a gauge theory 
that implements it ll22ll . The conserved scalar field suggests a 
generalization of electrodynamics. 

There is a sizeable literature on electrodynamics ll23i l24l 
25, 26] and non-abelian gauge theories |27, 2^|2pl on regu- 
lar lattices, mainly motivated by discretizing continuum theo- 
ries to facilitate numerical calculations. Chew ll2411 introduces 
a discrete vector calculus to formulate electrodynamics on a 
lattice. A Lagrangian approach has been used to obtain a net- 
work ap prox imation for electrodynamics in certain electronic 
devices [30]. An electronic-network model for coupled linear 
chemical reactions has also been proposed [31]. In this model, 
concentrations and not probabilities are mapped onto charges 
and a magnetic field or gauge fields are not introduced. 

The present Letter can also be read as a generalization of 
classical electrodynamics to a network and is thus of interest 
for the fundamental understanding of gauge theories, indepen- 
dently of the motivation from rate equations. It is not obvi- 
ous how one should generalize Chew's discrete vector calcu- 
lus l^^^l to a network, since there are no natural definitions of 
forward and backward differences. The Maxwell equations 
in integral form are more easily generalized. For example, the 
electric flux through a "wall" that divides the network into two 
parts is a meaningful quantity. It turns out that due to the lack 
of a length scale the electric (magnetic) field is the same as the 
electric (magnetic) flux and there is no difference between the 
differential and integral forms of the Maxwell equations. 

Maxwell equations. — Our goal is to construct electric and 
magnetic fields that implement the conservation of probabil- 
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ity. If we view Pi as a charge density, it should provide the 
sources of the electric field. We define an antisymmetric field 
Eij on the edges. Then Gauss' Law should read 

Y.E,,^P.,. (4) 

The sum over all outgoing electric fields equals the enclosed 
charge. Note that Eji is generally nonzero also if transitions 
from \i) to |j) are forbidden so that Jji = 0. This is not 
surprising since in standard electrodynamics there is certainly 
an electric field in an insulator As far as the electric (and, 
as we will see, the magnetic) field is concerned, the network 
is maximally connected, every node is the neighbor of every 
other node. 

Equation (|4]i leads to a problem: Summing over i we find 

i'j i 

We can rectify this by adding an additional fictitious node 
with Po = — 1 and RiQ = i?ot = for all i. Node does 
not affect the rate equations for the real states but makes the 
network charge neutral. The problem is that in Eq. (|5]) we 
obtain a relation between the flux through the surface of the 
system and its total charge. However, our network does not 
have a surface. The same problem arises if one tries to con- 
struct electrodynamics on any compact space. 

Next, we define a magnetic field B(ijk) on the elementary 
oriented plaquettes of the network, which are oriented trian- 
gles (ijk). The magnetic field is invariant under cyclic com- 
mutation of i,j,k and changes sign for anti-cyclic commuta- 
tions. To have Jy generate the magnetic field, we write the 
Ampere-Maxwell Law as 

1 . 1 

k c c 

Figure [Ha) shows that the equation contains a sum over the 
magnetic field penetrating all plaquettes adjacent to the edge 
carrying the current, corresponding to a line integral. Eji and 
B(^ijk) have the same dimension. Then c is a frequency, not a 
velocity, since the network does not have a length scale. 

Faraday's Law should relate the time derivative of the mag- 
netic field through a plaquette to the sum over electric fields 
along its edges. This is achieved by 

Eji + Ekj + Eik + - Bf^ijf.^ = 0. (7) 

The sum Eji + Ekj + Eik is the analogue of the electromotive 
force around the plaquette so that its inductance is 1/c. 

Finally, the magnetic flux through the surface of any vol- 
ume should vanish. An elementary cell on the network can be 
characterized by four nodes k, I. Then, 

B{tjk) + B(^jik) + B(^iik) + -S(jy) = 0. (8) 
Figure[Tlb) shows that the orientation of all faces is the same. 
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FIG. 1 : (a) Sketch showing that the Ampere-Maxwell Law relates the 
current on the edge from i to j to the circulation of the magnetic field 
around this edge, (b) Sketch of an elementary cell on the network, 
bounded by four plaquettes oriented outwards. 

The Maxwell equations imply the continuity equation: 
Summation of Eq. (|6]) over j gives 

jk j j 

The first term vanishes by symmetry. Using Eq. (|4|l we reob- 
tain Eq. Note that this derivation also works if the rates 
Rij depend on time. In fact, the derivation does not make any 
assumption about the dependence of Jij on Pi. 

Coulomb and Biot-Savart Laws. — We now consider the so- 
lution of the Maxwell equations for the fields. Direct calcu- 
lation yields P, - Pj = NEji + c"! J2k B(ijk) = NEj, + 
c^'^Eji + c^^Jji. Here, N is the number of nodes, including 
node 0. In the static case, we thus find the Coulomb Law 

E,^^^{P^~P,). (10) 

Interestingly, it is local: The static electric field on an edge is 
completely determined by the charges on the adjacent nodes. 

We also find c"i(J,-i + Jfej + Jjfc) = NBf^ijk) + c^^Bf^ijk). 
In the static limit we obtain a local Biot-Savart Law 

B{i-ik) = {J]i + -hj + Jik)- (11) 

Gauge fields. — As in standard electrodynamics, we intro- 
duce gauge fields to satisfy the homogeneous Maxwell equa- 
tions d?) and dHJ. The magnetic field is written as 

E{ijk) = Aji + Akj + Aik, (12) 

where Aij is antisymmetric. The form on the right-hand side 
is analogous to a curl, compare Eq. O. Interestingly, the mag- 
netic field has many more components than the "vector poten- 
tial" Aji, which is defined on edges, not plaquettes. 

Then Faraday's Law takes on the form Eji+Aji/c+Ekj + 
Akj/c + Eik + Aik/c — 0, which suggests to write 

Ej, = -{(f>j -<l>,)-^Aj,, (13) 

where the "scalar potential" (j)i is defined on the nodes. It is 
easy to show that Eqs. (|7]i and ([8]) are indeed satisfied. 
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The electric and magnetic fields are invariant under the si- 
multaneous gauge transformations 



A 



Aj^ + A-i - A,, 

Ai, 

c 



(14) 
(15) 



where A^ are arbitrary time-dependent functions. One can 
now consider specific gauge choices. For example, the ana- 
logue of the Lorenz gauge requires 



+ 0, = 0, 



0. 



(16) 



It is easy to show that this choice is possible. Then the re- 
maining two Maxwell equations take on the simple form 



P., 



1 



1 



NAji H — Aji — — Jji . 



(17) 



These equations are coupled by the gauge condition ( fTSI l. 

The representation in terms of gauge fields facilitates the 
discussion of the number of dynamic degrees of freedom. For 
a regular lattice. He and Teixeira lE6ll show that the electric 
and magnetic fields contain the same number of independent 
degrees of freedom. For the network, we first reformulate the 
question by asking how many electric and magnetic field com- 
ponents can consistently and independently be specified. 

We choose the gauge (pi = so that the Aij (t) are N{N — 
l)/2 independent functions. Since the Maxwell equations 
lead to second-order differential equations for the Aij (t), we 
can and must specify two initial conditions for each to deter- 
mine the solution. This gives N{N — 1) independent degrees 
of freedom. If we specify the initial values for the electric 
and magnetic fields instead, we must specify all N{N — l)/2 
components of Eij to fix Aij. We then must also specify 
N{N — l)/2 components of Bf^iji^^ to determine the solution. 
Thus the number of dynamic degrees of freedom in the elec- 
tric and magnetic field is the same also on the network. 

Lagrangian and energy. — The Maxwell equations can be 
concisely expresses by the Lagrangian 



L=l^El-l^ B, 



2 ^ 



2 

(ijk) 



E 



(18) 

As usual, Eij and -B(y7c) are to be expressed in terms of 
the gauge fields. The sum J2{ij) o^^r edges, counting 
each edge once, and J2{ijk) elementary plaquettes, 

counting each plaquette once. The action is iS = J dt L. It 
is straightforward to show that the Euler-Lagrange equations 
for Hamilton's principle SS = Q reproduce the two inhomo- 
geneous Maxwell equations. Incidentally, a co variant formu- 
lation is not possible due to the different structure of space 
(network) and time (continuum). 

The Lagrangian suggests to define the energy densities of 
the electric and magnetic fields as wfi = £'?,/2 and wf^^h = 
^(ijk)/'^' respectively. The local energy balance can then be 



expressed using the "Poynting vector" Sj^ 
With this definition we obtain 



-cEjiB^^jky 



d ci 
dt^'r^ 



di^iijk) - ^Ji ^kj Slf.. 



(19) 
(20) 



Energy is not conserved but changes due to ohmic dissipation. 
However, unlike in standard electrodynamics, the energy can 
actually increase, since Eji Jji can be negative for asymmetric 
rates Rij. Note that the generalization of the cross product 
yields a peculiar object S^^, which is symmetric only in its 
two lower indices. 

Medium equation. — The Maxwell equations do not yet 
form a closed set. As in standard electrodynamics in media 
one has to complement the Maxwell equations by medium 
equations describing the response of the medium. 

The situation most closely resembling our case is that of 
a conductor. Ohm's Law for the network would read — 
aij Eij. However, in our case the medium equation is just the 
definition of the current, Eq. (O. This equation is conceptu- 
ally different from Ohm's Law in that it expresses the current 
density in terms of the charges, not the electric field. The ori- 
gin is that we are actually describing a diffusive system. 

We can rewrite the medium equation to resemble Ohm's 
Law: With Gauss' Law ^ we obtain 



Jji — 2_^iPjiP'ki + RijEjk). 

k 



(21) 



This is a nonlocal version of Ohm's Law on the network, 
where the rates Rij play the role of conductivities. However, 
since Rij is not symmetric, nonzero currents can be present in 
the stationary state. 

If the rates are symmetric, Rij = Rji, which is a sufficient 
but not necessary condition for detailed balance lEoll . we can 
rewrite the medium equation as 



J 7' 



NRu E 



-E^ 

cN ^ 



(22) 



Here, NRji acts as the conductivity and the expression in 
parentheses is the force per unit charge. The first term is the 
normal electric field, while the second is not present in stan- 
dard electrodynamics. It is the rate of change of magnetic 
circulation around the edge from j to i. 

Discussion and conclusions. — We have formulated classi- 
cal electrodynamics on a network of possible states, motivated 
by the conservation of probability in systems of coupled rate 
equations. We close with a number of remarks. 

(i) Even if the transition rates are nonzero only between cer- 
tain states i, j, we have to introduce electromagnetic fields on 
all possible edges (ij) and plaquettes (ijk) of the network. 
Thus this system does not reduce to electrodynamics on a lat- 
tice 1 24] for the case that the nodes connected by allowed 
transitions form a regular lattice. Nevertheless, the number 
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of degrees of freedom of the electric and magnetic fields on 
the network are equal, as for a regular lattice. 

(ii) As in standard electrodynamics, only the two inhomo- 
geneous Maxwell equations are required to implement the 
continuity equation. Thus the theory allows to introduce mag- 
netic monopoles, which here live on the dual network formed 
by the cells of the original network. 

(iii) Compared to continuum electrodynamics, scalar fields 
correspond to quantities defined at the nodes (Pi, (t)i), polar- 
vector fields correspond to quantities defined on edges (J^ , 
Eij, Aij), and axial-vector fields correspond to plaquette 
fields (Bj-yTj-j). A Poynting vector can be defined, but has 
a more complicated structure than the "vector" Eij because 
electric and magnetic fields live in different places. 

(iv) We have reformulated the rate equations as a variational 
principle, SS = 0, for the action of electromagnetic fields. 
What is gained by this formulation? On the one hand, it rep- 
resents a new way to think about rate equations within the 
framework of electrodynamics. It works for any dependence 
of the current on the probabilities and their time derivatives, 
as long as Jy is antisymmetric. 

The framework is expected to be useful for inverse prob- 
lems: How can one construct a dynamical system with spec- 
ified properties? As an example, note that the new fields 
E'j, = Ej, + j:kB(^jk)/cN, B[^^,^ = and the new 

cun-ent Jj, = J^i - cY.kB(iok) - T.kB(i]k)/cN satisfy 
the Maxwell equations with the original charges Pi. Since 
B' = 0, Jjj + + J-j, vanishes for all i, j, fc, cf. the dis- 
cussion leading to Eq. ( fTTI ). This means that for any time- 
dependent probabilities Pi one can find a dynamical system 
with these probabilities and with Jji + Jkj + Jik — 0. Since 
the current is curl-less, one can write it in terms of a "current 
potential," Jji = ~{tpj — tpi). Gauss' Law then shows that 
4>i = —Pi/N is a solution. 

(v) The framework is of interest beyond the motivation in 
terms of rate equations, since it shows that a classical gauge 
theory can be formulated consistently on a network, which 
does not have a natural metric. For this reason the parame- 
ter c is a frequency instead of a velocity. It is an interesting 
question how non-abelian gauge theories fare in this context. 

The author would like to thank R. K. P. Zia, J. M. Arnold, 
F. S. Nogueira, and J. P. Ralston for illuminating discussions. 
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